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ABSTRACT 

An equation of state for the domain extending from hot gases to cool-dense fluids 
is formulated for a hydrogen-helium mixture. The physical processes take account of 
temperature ionization and dissociation, electron degeneracy, Coulomb coupling and 
pressure ionization. Pressure ionization and Coulomb coupling are studied with sim- 
ple and comprehensive modeling. A single and complete algorithm is achieved with 
explicit expressions available for the whole domain from hot gases to cool dense fluids 
(T > 10 2 K). Pressure ionization and Coulomb coupling have been examined for their 
contributions to the pressure and internal energy. The result reveals that their contribu- 
tions smooth the variation of the pressure and internal energy in the region of pressure 
ionization even at very low temperatures. 



Subject headings: atomic processes, equation of state, stars: low mass, planets 



1. Introduction 

The equation of state (EOS) of cool dense fluids is concerned with the interiors of very low mass 
stars, brown dwarfs and giant planets in astrophysics. The complexity of the physical processes of 
such an EOS stems from the inclusion of interparticle interactions among atoms, molecules, ions 
and electrons. The physical processes concerned mainly include two types, i.e. pressure ionization 
and Coulomb coupling. Many studies have been carried on for Coulomb coupling among charged 
particle. For instance, the internal energy of fully ionized plasmas of classical one-component 
has been given by advanced physics models (see Slattery and Doolen 1980) for the case of strong 
coupling. The free energy of fully ionized hydrogen-helium mixture has also been evaluated for 
thermodynamic properties with strong coupling (see DeWitt and Hubbard 1976). The Coulomb 
coupling in weak limits can be well described by Debye-Huckel model (see Cox and Giuli 1968; 
Mihalas et al. 1988). However, the Coulomb coupling between strong and weak limits still remains 
for the study of its model and effects. In addition, the problem of pressure ionization has been 
remaining thorny in the EOS calculation for cool dense fluids (see Saumon et al. 1995). Although 
the advanced physical models of pressure ionization have been studied for the EOS (see Blenski 
and Ishikawa 1995; Saumon and Chabrier 1991, 1992), it is not easy to incorporate them into the 
astrophysical EOS for their practical use and investigation. In this paper, we establish the physical 
models for pressure ionization and Coulomb coupling which have simple expressions and can be 
easily incorporated into the EOS formalism. 
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Almost all of the EOS for cool dense fluids are given with their results in tabular form for 
astrophysical application (e.g. Fontain et al. 1977; Magni and Mazzitelli 1979; Mihalas et al. 
1988; Saumon et al. 1995). This is because the complex formalism and time-consuming calculation 
of the EOS for cool dense fluids stem from many physical processes. They include temperature 
dissociation and ionization, electron degeneracy, Coulomb coupling and pressure ionization, etc. 
However, explicit and simplified expressions for a comprehensive EOS can be of benefit to the study 
of the unfamiliar physical processes concerned in the EOS and to the study of their importance 
to and influences on astrophysical phenomena. We aim to present such an algorithm for the EOS 
covering the domain from cool dense fluids to hot gases. 

Furthermore, it is also important to assure the reliability for the formulated EOS and its 
results. Therefore, we not only establish the explicit and simple expressions for the EOS but also 
compare the resulted EOS with those from other more sophisticated models. In order to preserve 
the continuity and generality of the EOS, we formulate the EOS on the basis of the grand canonical 
approach in the chemical picture for the EOS (see Luo 1994). It allows to establish a single algorithm 
of the EOS available for the whole domain from cool dense fluids to hot gases, even when many 
physical processes are taken into account. The grand canonical approach in the chemical picture 
also allows to include mixed elements in the EOS. The EOS presented in this paper is formulated 
for a hydrogen(H)-helium(He) mixture and takes into account the physical processes of atomic 
ionization, molecular dissociation, electron degeneracy, Coulomb coupling, pressure ionization, etc. 

In Sect. 2 we establish the models of pressure ionization and Coulomb coupling. Sect. 3 presents 
an algorithm of calculating thermal equilibrium. The algorithm of calculating the pressure and in- 
ternal energy is presented in Sect. 4. In Sect. 5, the calculated results of the EOS are compared with 
other EOS, and are examined for the contributions of pressure ionization and Coulomb coupling to 
thermodynamic functions. 



2. Pressure ionization and Coulomb coupling 

We start from a H-He mixture gas with mass density p, temperature T and the mass abun- 
dances X and Y for H and He. In a volume V, the number of nuclei of H and He are given 
by 

N n = N A pVX/A n , N Uc = N AP VY/A Hc , (1) 

with the Avagadro's number Na, the atomic weights Ah = 1.0079 and ^He = 4.0026. 

Two physical processes, i.e., Coulomb coupling of charged particles and pressure ionization of 
atoms, have been taken into account for the nonideal effects on the EOS. Owing to the complexity 
of such physical processes, their incorporation into the EOS requires some assumptions to allow 
simplified and approximate models. 



2.1. Coulomb coupling 

The Coulomb coupling has considerable influences on the EOS, especially in fully ionized 
plasmas. In partially ionized fluids, it plays a subtle role in the determination of the ionization 
degrees as well as thermodynamic quantities. We first approximate the Coulomb coupling by a 
one-component model of fully ionized plasmas, and then add the corrections of partial ionization 
and mixed elements to its formalism. 
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The Coulomb coupling can well be described in two limit cases, i.e., weak coupling (r <C 1) 
and strong coupling (T 3> 1). The Coulomb coupling parameter V = (Ze) 2 /fkT is denned by the 
nuclear charge Z, the ion sphere radius f and temperature T. For a H-He mixture, we define the 
average nuclear charge Z as 

= jVH + 2jV H e 

The ion sphere radius f is determined by the number of ions AT ion in the volume V using the 
approximation (47r/3)f 3 iV"i on /y = 1. Considering the correction from the partial ionization, we 
approximate the ion number by AW = N e /Z with the number of free electrons iV c . As a result, 
the Coulomb coupling parameter can be written as 

{Zef (AkN c \^ 



kT V 3Zy 

with the corrections from the partial ionization and mixed elements. 

In the case of strong coupling (r > 1), the energy of a fully ionized plasma with one component 
(of nuclear charge Z) is found to be 

E cc = -kTN ion x 0.89752r, (T » 1) (4) 

with the ion number Ni on (see Slattery et al. 1980). For weak coupling r < 1, the Debye-Hiickel 
model can be used to obtain its Coulomb coupling energy (see Cox & Giuli 1968). We ignore the 
effect of electron degeneracy and suppose a fully ionized plasma with one component (of nuclear 
charge Z). We thus have 

E cc = — kTN[ on x (3/4) 1 / 2 [(1 + l/Z) vf 2 , (r « 1). (5) 

The results for both the limit cases indicate that the Coulomb coupling energy per ion in unit kT 
(i.e. E cc /N ion kT) is a function of the Coulomb coupling parameter V and the nuclei charge Z. The 
functions (E cc /N ion kT) of the two limit cases for Z = 1 are plotted versus T in Fig. 1 which shows 
their intersection at T = 0.13. 

In our EOS formalism the average potential energy of each free electron from Coulomb coupling 
e cc is adopted. The internal energy from the Coulomb coupling is calculated in terms of 

which is the result of the EOS with the grand canonical approach in the chemical picture (see Luo 
1994, 1996). With this definition and the constraints of the results of Coulomb coupling energy in 
two limit cases, the average energy of each free electron from Coulomb coupling can be fitted to be 

e cc = -kT 9 -^! (7) 

with a fitting function g(F, Z) defined by 

l/g = 1/0! + 1/02, (8) 
9l = (2/3)(3/4) 1 / 2 (l + l/Z) 3 / 2 r 3 / 2 , 
g 2 = 0.89752r . 
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The internal energy from Coulomb coupling resulted from the fitting function is therefore given 

by 

aTnr) • (9) 

according to Eqs. (3), (6) and (7). It has asymptotic form as in equation (4) for strong coupling 
case and as equation (5) for weak coupling case. The results from equation (9) for the Coulomb 
coupling energy per ion in unit kT {E cc /N\ ori kT) with Z = 1 are plotted in Fig. 1. It is compared 
with the results of the two limit cases. The accurate results for T > 0.3 (see Slattery and Doolen 
1980) are also plotted in Fig. 1. They indicate that the simply fitted formula still overestimate the 
Coulomb coupling energy about tens percent in the intermediate range of T. However, it is easy to 
add the power weighting to the fitting formulae to obtain better fitted results. 



Fig. 1. — The fitting function for logarithmic (E/kTNi on ) versus logarithmic T (solid line) , with those of 
the weak Coloumb coulping T <C 1 (dotted line) and of the strong Coloumb coupling r> 1 (dashed line) , 
and with the accurate results for T > 0.3 (circles) given by Slattery and Doolen (1980). 

The fitted formula of equation (7) has a derivative with respect to lnT 

d^g = 5 _ (3/2)ff 1+ g 2 

dlnT 2 gi +g 2 ' 1 ' 

with the functions g\ and g 2 defined in equation (8). It has asymptotic values 3/2 for V <C 1 and 1 
for T S> 1. It will be used in Sect. 4 for calculating the pressure and internal energy from Coulomb 
coupling. 



2.2. Pressure ionization 

Pressure ionization is caused by the interaction of an atom with surrounding particles which 
is responsible for the dissolution of atomic energy states. In the gases with high temperatures (e.g. 
T > 10 5 K), a great fraction of free electrons and ions lead to the Stark ionization of atoms. This has 
been well described by the MHD EOS with the occupation probability formalism (see Hummer & 
Mihalas 1988). For cool dense fluids, however, neutral atoms and molecules dominate the interaction 
on an acted-upon atom. Accordingly, the sophisticated physical models are established to produce 
continuous solutions of the EOS (e.g. Saumon & Chabrier 1992; Blenski & Ishikava 1995). These 
models are based on the advanced physics theory and partly examined by experimental data. 
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For astrophysical application and research, however, we may make approximate and statistical 
models for the complicated physical processes of pressure ionization. This will allow to examine 
the uncertainty of pressure ionization and its influence on the models of astrophysical objects 
through the EOS. 

Our model for pressure ionization is based on the energy change of an acted atom by the 
pressure impact of acting particles. When an acted-upon atom is collided by other particles, the 
pressure impact can always increase its energy. If the acted-upon atom is supposed to have an 
effective volume vi, the impact by acting particles with an equivalent pressure P\ will increase the 
energy of the acted-upon atom by an amount of AE\ = P\V\. This is the result from the view of 
statistical mechanics. In thermodynamic equilibrium with temperature T, the energy change AEi 
leads to a Boltzmann factor W = exp(— AEi/kT) for the existence probability of the atom. We 
thus can use the resulted factor W to present the occupation probability for the specified atom. 

The acting particles include all classical particle with the number of atomic neclei Nn + Nn e . 
They also include the bound electrons of other atoms beyond the mean separation d that the 
acted-upon atom can feel. The mean separation of atoms d is defined by 

4tt 3 fN H + N Ue 



3 d W^J =1 < u ) 

for a H-He mixture with the number of atomic nuclei Nu + Nn e . 

Suppose the atoms of H and He have the hydrogen-like wave functions of ground states w(Zr) = 
\^2Zre~ 2Zr , we then have the number of the bound electrons beyond the mean separation d 

iV c b = N H / w ( r )dr + 2Ar Hc / w(2r)dr = N H f(d) + 2N Be f(2d), (12) 

Jd Jd 

with a function f{x) defined by 

f(x) =e- 2x (l + 2x + 2x 2 ). (13) 

It is a rough approximation that the mean separation of atoms is used to determine the number 
of the acting electrons. Because an atom with a kinetic energy (3/2)kT can always approach another 
atom closer than the mean separation d, it can feel more bound electrons of other atoms than those 
outside the separation d. As a result, the integration limit d in equation (12) should be corrected by 
taking account of the thermal motion of the atom. Nevertheless, the thermal motion of the atom 
play a role in the determination of the integration limit at relatively high temperatures. In that 
case, temperature ionization play a more important role than pressure ionization. We find that 
the inclusion of the thermal motion dose not have a very important effect on the model of pressure 
ionization in the EOS formalism. In addition, its inclusion will complicate the model of pressure 
ionization in the EOS formalism. We therefore do not present its formulation in this paper. 

The acting atoms and ions are classical particles while the acting electrons with the number 
behave like the Fermi-Dirac particles. As a result, the equivalent pressure P\ from the acting 
particles is given by the Fermi-Dirac statistics 

/ AT H + A/wp \ 2 N h F V2 (\i) 

with the Fermi-Dirac integral F 3 / 2 and Fy 2 - The degeneracy of the acting electrons Ai is determined 

N*/V = C e T*' 2 F 1/2 (\ 1 ) (15) 
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with a constant C e = (y/2/7r 2 )(km c /Ti 2 ) 3 / 2 (see Luo 1994). 

The effective volume v\ of the acted-upon atom is defined by its radius rj in terms of v\ = 
(47r/3)r|. As a result, we obtain the occupation probability in terms of W = exp(— P\V\/kT) in a 
logarithmic form 

(16) 

for an atom with an effective radius r y In equation (16) the first term is identical to the result of 
the hard sphere model of atoms (see Hummer and Mihalas 1988). The second term is the result 
of the interaction from the bound electrons of other atoms on the acted-upon atom. The radii of 
atomic hydrogen (HI) and singly ionized helium (Hell) are taken to be rm = ao an d rneii = 0.5ao 
with the Bohr radius ao, and that of atomic He is assumed to be rn e i = 0.8ao- The resulted 
occupation probabilities are thereby represented by Whi, Wud and Hjieii- 

The occupation probability given by equation (16) is an explicit function of T, Ah and A^e as 
well as Tj. It thus can be incorporated into the calculation of the thermal equilibrium. Its simplified 
form allows to obtain analytical expressions for their contribution from pressure ionization to the 
pressure and internal energy (see Sect. 4.5). 



4-7T q 

lnW f . = __ r 3 



3. Thermal equilibrium quantities of the EOS 

3.1. Thermal equilibrium quantities 

Thermal equilibrium quantities are explicitly defined by the degrees of atom ionization and 
molecule dissociation. For the EOS in the chemical picture, the determination of thermodynamic 
functions significantly depends on the degrees of ionization and dissociation. According to the 
grand canonical approach in the chemical picture (see Luo 1994), the thermal equilibrium of a 
mixture can be specified by the probabilities of atomic nuclei in various configurations. 

We total include six configurations for H and He. They are specified for hydrogen by molecular 
state H2 (abbreviated as M), atomic state HI and ionized state HII, as well as for helium by atomic 
state Hel, singly ionized state Hell and doubly ionized state Helll. The probabilities of hydrogen 
nuclei in the configurations H2, HI and HII are thereby specified by xm, xm and xhii, with a 
normalization condition xm + xhi + xhii = 1- The probabilities of helium nuclei in the configurations 
Hel, Hell and Helll are specified by xh g i, XHell an d XHell with xn e i + XHcii + XHciii = 1- 

The energy of an atomic configuration is equal to the energy of bound electrons in it, and 
is determined by the ionization energy and dissociation energy. For hydrogen and helium, their 
configuration energies are defined by 

e h 2 = -2xm - XH 2 
Em = -Xhi 
Emi = 

Eftel = — XHel ~ XHell 
-^Hell = — XHell 

^Helll = (17) 

The energies of ionization and dissociation xh 2) Xhi, XHel, XHell for H and He have values 4.48, 
13.598, 24.587, and 54.416 in eV. The energy of a fully ionized configuration is equal to zero because 



- 7- 



it has no bound electron. The configurational energy of H2 specified by Eu 2 is equal to the energy 
of two electrons bound in H2. 

The probabilities of atomic nuclei in various configuration are alternative quantities for de- 
scribing ionization and dissociation degrees. The number of atomic nuclei in each configuration is 
equal to the number of atomic nuclei multiplied by the corresponding probability. The number of 
free electrons is thus given for a H-He mixture by 



A H-He mixture in thermal equilibrium is described by its temperature T, volume V and the 
number of hydrogen nuclei Ah and that of helium nuclei An e - Following the conservation equation 
for the electron number (see Luo 1994) we have an equation for thermal equilibrium 



Of equation (19) the left-hand side counts the total number of electrons from hydrogen and 
helium while the right hand side represents the distribution of electrons in molecular, atomic, ionic 
and free states. On the right-hand side, the first term is the number of free electrons, the second 
term counts the bound electrons in molecular and atomic configurations of hydrogen, and the third 
term counts those in the atomic and ionic configurations of helium. 

Equation (19) has an unique unknown variable which is the electron degeneracy A. The 
probabilities of nuclei in various configurations xm, xhi, XHci and XHeii are all its explicit functions 
according to Sects. 3.3 and 3.4. The Coulomb coupling energy for each free electron e cc is defined 
in Sect. 2.2. It is also an explicit function of the electron degeneracy A by means of the free electron 
number N e . As a result, the equation for thermal equilibrium, equation (19), can be solved for the 
electron degeneracy A as well as other equilibrium quantities xm, £hi, xh c i, £Hcii and N e . 



Helium has two degrees of ionization and thus has three configurations Hel, Hell and Helll. 
The probabilities of atomic nuclei in various configurations can be determined by the formalism of 
the atomic configuration probability (see Luo 1994). The grand partition function of electrons in an 
atomic configuration \NE > is defined by Z = gW exp(AA — E/kT) with the statistical weight g, 
the number of bound electrons A, the configuration energy E (i.e. the energy of bound electrons) 
and the occupation probability W. 

As a result, we have the grand partition functions for various configurations of helium 



A c = A H x H ii + A H c(xhcII + 2x Hc iii) • 



(18) 



3.2. The equation for thermal equilibrium 




3.3. Equilibrium quantities of helium 



Z-Hel = V 
Zu&l = 2 
^HcIII = 1 



W ne i exp[2A - E Hel /kT] 
2W He iiexp[A- E HeI1 / 'kT] 



(20) 
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which result from g = 1 and N = 2 for Hel and g = 2 and N = 1 for Hell. The fully ionized 
configuration Helll has its grand partition function equal to unity because it has no bound electron 
and is characterized byA r = 0,i? = 0,(/ = l and WHeiii = 1- The occupation probabilities Whi 
and W^Heii which account for pressure ionization, are formulated in Sect. 2.2. 

The probabilities of helium nuclei in various configurations are therefore given by 

XHcI = ^Hel/^He 

XHell = ^Hell/^He 

XHelll = ^HcIIl/^He (21) 

2-Re = ZucI + -Zfiell + ^Helll 

where Zu e is the total grand partition function of elections for He. They satisfy a normalization 
condition x He i + XHell + x Hc iii = 1. 



3.4. Equilibrium quantities of hydrogen 

We include three configurations H2, HI and HII for hydrogen in the EOS formalism. Because 
of the inclusion of molecular hydrogen the formalism of the probability of hydrogen nuclei has a 
different form from those of helium. We first give the probabilities of hydrogen nuclei in HI and 
HII for which no molecular state is taken into account 

ym = 2>m/ 

ymi = Zmi/Zn (22) 
Zyi = Zyii + Zmi 

with normalization ym + ynn = 1- The grand partition functions of electrons in HI and HII are 
defined by 

Z m = 2W m exp[\-E m /kT] 

Zmi = 1 (23) 
which are the results of g = 2 and N = 1 for HI, and the fully ionized configuration for HII. 

When we introduce the probability of hydrogen nuclei in the molecular configuration xm , the 
exclusion of hydrogen nuclei in the molecular configurations leads to the actual probabilities of 
hydrogen nuclei in configurations HI and HII 

xm = ym(l - xm) 

xmi = ymi(l - (24) 

For the dissociation of molecular hydrogen, the pressure equilibrium constant defined by Vardya 
(1960) leads to 

K (T) = ?t= ( fcTiV Hi/^) 2 = 4i 2fcriY H 

1 ; P H2 kTN n jV x M V { ' 

with the probabilities of hydrogen nuclei xm m H2 and xhi in HI. Equation (25) has employed the 
number of hydrogen molecules A^h 2 = -^hxm/2 and the number of hydrogen atoms A^hi = -/Vhxhi- 
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With the normalization conditions xm + %m + ^hii = 1 and ym + Uhu = 1, Eqs. (24) and (25) 
lead to a solution for the probability of hydrogen nuclei in the molecular configuration 

x M = 1 + R - VR 2 + 2R (26) 

with 



^^(TVh/F)' 

The Vardya's (1960) expression for the pressure equilibrium constant can be written as 

K{T) = P exp[-(Z>i/fcT) + (D 2 /kT) 2 - (D 3 /kT) 3 } (28) 

with a pressure constant Po = 3.4159 x 10 12 dyn/cm 2 . The constants Di, D 2 and D3 take values 
4.92516, 0.237047 and 0.148409 in eV. 

Temperature ionization and pressure ionization of hydrogen atoms have fatal effects on the 
dissociation of molecular hydrogen. The reduction of the number of hydrogen atoms rapidly de- 
creases the number of hydrogen molecules. In the above formalism, the decrease of yni reduces xm 
in terms of xm ~ 1/2R for R S> 1 according to Eqs. (26) and (27). 



3.5. Thermal equilibrium calculation 

We finally established an algorithm for calculating thermal equilibrium quantities from Eqs. 
(19)-(28). Equations (21), (24) and (26) define the probabilities of atomic nuclei in various con- 
figurations as functions of the degeneracy A. The equation for thermal equilibrium, equation (19), 
therefore becomes an enclosed equation with one variable A for which can be solved. The iteration 
method is an effective and fast way to solve equation (19) in most region of the p — T domain 
except in the region where strong pressure ionization takes place. In contrast, the bisection method 
is powerful and can always arrive at unique convergent solutions in the whole p — T domain. 

The thermal equilibrium of a mixture can be well revealed by its mean molecular weights in a 
p — T domain. The mean molecular weight p is defined as the mean weight per particle in atomic 
mass unit, i.e. p = pVN\/N with the total number of particles N in the mixture. We thus have 
the mean molecular weights p for a H-He mixture in terms of 

IX Y 

- = -j-(l - X M /2 + Xffll) + — (1 + XHcII + 2x He iii). (29) 
P Aft Aft e 

It is the result of the degrees of ionization and dissociation through the probabilities of atomic 
nuclei in various configurations. 

We performed a calculation of thermal equilibrium for a H-He mixture with abundances X/Y = 
0.7/0.3. The reciprocals of the resulted mean molecular weights 1/p are plotted in Fig. 2 for a p—T 
domain. It shows three slopes for low densities log p < from low temperatures to high ones. They 
correspond to the dissociation of hydrogen molecules, the ionization of hydrogen atoms and that 
of helium atoms respectively. There are other three steep slopes in the density range < log p < 2 
for low temperatures logT < 5, which are the results of pressure ionization of hydrogen atoms and 
helium atoms. The surface plot of Fig. 2 shows that the minimum and maximum values of 1/p, 
0.422 and 1.614, correspond respectively to the gases of completely neutral helium and molecular 
hydrogen and those of completely ionized mixture. 
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Fig. 2. — The reciprocals of the mean molecular weight 1/p in a p — T domain for a H-Hc mixture 
X/Y = 0.7/0.3. 

In addition, Fig. 2 reveals that the densities for pressure ionization of hydrogen are located 
between 1 to 10 g/cm 3 at very low temperatures while those of helium between 10 to 100 g/cm 3 . The 
densities for pressure ionization are almost independent of temperatures, and pressure ionization 
take place in very narrow density ranges. However, we realize that the occupation probability 
by equation (16) is formulated especially for cool-dense fluids. Although it can produce pressure 
ionization for partial ionized gases at high temperatures, it underestimates the effect of pressure 
ionization by the Stark effect for logT > 4.5. 



The EOS uniquely determines the pressure P(p, T, X) and internal energy E(p, T, X) as well as 
other thermodynamic quantities as functions of density p, temperature T and element abundance 
X, Y. The pressure P, the internal energy density per unit mass E/pV and their derivatives with 
respect to p and T are required to dictate the mechanical and thermal equilibria of stars and 
planets. In this section, we present the comprehensive expressions for the pressure P and internal 
energy E as functions of T, Nr and A^e which are related with T, p and X by equation (1). 

We include the following contributions to the pressure and internal energy: (1) the classical 
particle (molecules, atoms and ions) and the free electrons as Fermi-Dirac particles; (2) Coulomb 
coupling for P cc and E cc ; (3) pressure ionization for P p j and E p [. The internal energy also include (4) 
the configurational terms of atoms, ions and molecules E con f and (5) the vibrational and rotational 
terms of hydrogen molecules E TV . 

As a result, the pressure and internal energy for a H-He mixture is written as 



4. Thermodynamic functions of the EOS 



4.1. Pressure and internal energy 




"I" Pec ~\~ Ppi 



(30) 



E = 




+ E cc + E pi + -E con f + E, 



rv 



(31) 
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with the number of the classical particles N c = Nn(l — XM./2)+Nn e and the number of free electrons 
-/V e given by equation (18). In the Fermi-Dirac integrals F\i<i and -F3/2, the effective term of the 
electron degeneracy A e = A — e cc /kT is determined by the electron degeneracy A and the Coulomb 
coupling energy for the free electron e cc . The first and second terms in Eqs. (30) and (31) represent 
the ideal contributions from the classical particles and the free electrons. 



4.2. Energy of configurations 

The configurational energy comes from the energy of bound electrons in molecules, atoms and 
ions. This energy is determined by the numbers of nuclei in various configurations and the energies 
of various configurations are defined by equation (17). For a H-He mixture, the total energy from 
various configurations is given by 

-Econf = N h (x m Eh 2 /2 + xmEm + xmiEmi- E H2 /2) 

+ ^Hc^Hel-ElM + XHell-^Hell + ^Helll-EHelll — -Erci) (32) 

The last terms in two brackets, i.e. (— En 2 /2) and (— End), are added to retain a zero point of the 
configurational energy in the case of 100% molecular hydrogen and neutral helium. 



4.3. Rotational and vibrational energies of hydrogen molecules 

The rotational and vibrational states of hydrogen molecules (H2) have contribution to the 
internal energy. The characteristic temperature of its rotation is r = 85.38K and that of its 
vibration is @ v = 5987K (see Gopal 1974). For temperatures T> 6 r which is our interest for cool 
dense fluids, the rotation of H2 is sufficiently excited and its rotational energy reaches the classical 
value of kT. In contrast, the vibration of hydrogen molecules cannot be excited to classical states 
at temperatures T < 6 V or T ~ 6 V . Therefore, we use the Einstein function (see Gopal 1974) for 
the vibrational energy of H2 . 

As a result, the contribution from the rotation and vibration of H2 to the internal energy is 
given by 

t9 » +1*6,1 



E TV = kTN R2 + 



N R2 (33) 



_ e e v /T_ 1 2 

with the number of hydrogen molecules Nn 2 = Nhxm/2- In equation (33) the first term results 
from the rotation of H2 while the second term from its vibration. The factor E IV /kTNn 2 is sensitive 
to temperatures around T ~ V . 



4.4. Contribution from Coulomb coupling 

The contribution from Coulomb coupling to the pressure and internal energy is calculated from 
the Coulomb coupling energy contributed to the free electrons. Equation (7) defines the average 
energy from Coulomb coupling for each free electron e cc . According to the formalism of the EOS 
formulated by Luo (1994), the pressure from Coulomb coupling is given by 

_iVe / de cc \ __ iVeecc 1 ( d\ng \ 
cc ~ V\dlnVj TN ~ V 3UlnI7 1 j 
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and the internal energy by 

^-^-"•G^L-MIe?) (35> 

which are the results of equation (7). The derivative factor in Eqs. (34) and (35) has been given 
by equation (10) and is a function of the Coulomb coupling parameter T and the average nuclear 
charge Z. A relation P cc = E CC /3V is retained for the Coulomb coupling from Eqs. (34) and (35). 



4.5. Contribution from pressure ionization 

Pressure ionization is caused by interparticle interactions and has its contributions to the 
pressure P and internal energy E. Its contributions have been formulated as the results of the 
occupation probability W of bound atomic configurations (see Luo 1994) in terms of 



1 



Ni fdlnWi 



V V dhxV 



(36) 



E, 



pi 



2 4^ 3 V dlnT 



(37) 



VNi 



where the sums go over all bound configurations for j = {H2, HI, Hel, Hell} with the species number 
Nj. The factor (1/2) is introduced to take account of the fact that the interparticle interactions 
responsible for pressure ionization have actually been counted twice for each atom in the summation. 

The fully ionized ions, e.g. HII and Helll, have no contribution because they have no bound 
systems interacting with surroundings. The contribution from molecular hydrogen should be in- 
cluded in that from hydrogen atoms since its pressure dissociation results from pressure ionization 
of atomic hydrogen according to Eqs. (26) and (27). As a result, the contribution from pressure 
ionization to the pressure for a H-He mixture is given by 



pi 



+ 



V 

NHe 

V 



dlnW m \ 



XHel 



dln^ 
dlnW Hel 
<91n^ 



+ a^Heii 



'dlnWkcn 
v dlnV , 



(38) 



and that to the internal energy given by 



E pi = -kTN H (x M + xm)( 



+ 



-kTNn c 



XHcl 



d\nW m \ 
d\nT ) 
d In V^Hei ' 



<91nT 



+ ^Hell 



dlnWucii 
<91nT 



(39) 



The derivative terms in Eqs. (38) and (39) are determined by the formalism of the occupation 
probability Wj. They have analytical expressions for j = {HI, Hel, Hell} 



d In Wj 
dlnV 



47T n 

TNj 3 3 



(40) 
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d In Wj\ 4vr 



<91nT 



—rp^F^X^ihiX,)-!] (41) 



as the results of equation (16). The degeneracy parameter of the acting electrons Ai is defined by 
equation (15). The factor h(Xi) is a function of Ai and has a form 

h A ) - 3i? V 2(Al) (42) 
^3/2(^1)^-1/2(^1) 

which has asymptotic values h = 1 for Ai <C and h = 5/3 for Ai 3> 0. 



4.6. Discussion 

In this section, we archive a complete algorithm for calculating the pressure and internal 
energy with Eqs. (30)-(42). All of the expressions have explicit and simple form although some 
approximations have been adopted. The first approximation is the use of one-component plasma 
model for Coulomb coupling in Sect. 2.1. However, it is a good approximation when the corrections 
from partial ionization and element mixture are added. And it is valid for a very wide range of the 
Coulomb coupling parameter with < F < 160. The second one is the model of pressure ionization 
established in Sect. 2.2. It is based on the concept of the pressure impact of acting particles on 
the acted-upon atom, and represents the ionization of atoms by the high pressure of matter. Its 
simplicity allows to retain the analytical expressions for the contribution from pressure ionization 
to the pressure and internal energy. 

The algorithm of the EOS established in Sect. 3 and 4 is based on the chemical picture for 
mixed elements. It allows to easily include more elements for their contribution to the EOS. The 
formalism for hydrogen can be developed for the elements with molecules while that for helium can 
be developed for the elements with multielectron atoms. 



5. Results 



5.1. Contribution from nonideal effects 

Pressure ionization and Coulomb coupling lead to the nonideal effects on the EOS. They not 
only increase the ionization degrees, but also contribute to thermodynamic functions. In Sect. 4, 
the formalism for the EOS separates their contribution terms from those ideal terms of classical 
particles and free electrons. This allows to examine the importance of pressure ionization and 
Coulomb coupling for their contribution to the pressure and internal energy. 

In spite of the reliability of our EOS for very cool-dense fluids, we performed the calculation 
upto very high densities for low temperatures (T > 10 2 K). That reaches the states of molecular 
solid and metallic solid (see Kerley 1972). The simplicity of the EOS formalism allows to easily 
obtain the solution and to roughly investigate the properties of matter from hot gases, to cool-dense 
fluids, and to solids. 

We first examine the contributions from pressure ionization and Coulomb coupling to the 
pressure and internal energy for a H-He mixture with X/Y = 0.7/0.3. Their fractions over total 
pressure and internal energy are plotted in Fig. 3 and Fig. 4 for a p — T domain. They are actually 
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formulated by (P pi + P CC )/P and (E pi + E cc )/E . Equations (34), (36), (38) and (39) indicate that 
pressure ionization has positive contributions to the pressure and internal energy while Coulomb 
coupling has negative contributions to them. 



Fig. 3. — The fractions of the contribution from pressure ionization and Coulomb coupling over the total 
pressure in a p — T domain for a mixture 



Fig. 4. — Same as Fig. 3, but over the total internal energy. 

Figures 3 and 4 reveal that the contribution from pressure ionization increases the pressure P 
and internal E before pressure ionization takes place. In contrast, the contribution from Coulomb 
coupling reduces P and E when pressure ionization occurs and a great fraction of charged particles 
are produced. In addition, the contribution from pressure ionization dominates both pressure and 
internal energy in a density range of — 1 < logp < 1 for temperatures logT < 4.5. The contribution 
from Coulomb coupling can reduce the pressure and internal energy by tens percent just when 
pressure ionization takes place at very low temperatures. 

In Figs. 5 and 6 the pressure P and internal energy density E/Vp are plotted on isotherms 
in a density range — 2 < logp < 2 in solid lines. Together plotted in dotted lines are the ideal 
terms, i.e., P — Ppi — Pec and (P — Ppi — Pec) P which are the pressure and internal energy 
excluding the contributions from pressure ionization and Coulomb coupling. From Fig. 5 and Fig. 
6 we see the following points. (1) The densities of pressure ionization are located in a very narrow 
range. There the ideal term of free electrons increases the pressure and internal energy rapidly and 
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significantly. This is the result of the fast production of free electrons due to pressure ionization. 
(2) The contribution from pressure ionization starts to play a role and increases the pressure and 
internal energy far before pressure ionization takes place. This results from the interaction among 
atoms when they are close to each other. (3) The contribution from Coulomb coupling plays a 
role just after pressure ionization takes place. A great fraction of charged particles are produced 
and the resulted Coulomb coupling from them reduces both pressure and internal energy. (4) The 
increase of the pressure and internal energy by the contribution from pressure ionization before 
occurring of pressure ionization and the decrease by that from Coulomb coupling after occurrence 
of pressure ionization smooth the variation of both pressure and internal energy in the region 
of pressure ionization, in spite of the fact that the fast production of free electrons by pressure 
ionization increase the ideal pressure and internal energy steeply and significantly. 



Fig. 5. — Pressure isotherms (in unit dyn/cm 2 ) in solid lines, compared with those of the ideal terms in 
dotted lines. The isotherms are (from bottom to top): logT = 2.0,2.5,3.0,3.5,4.0,4.5,5.0,5.5,6.0. 



Fig. 6. — Same as Fig 5 but for internal energy density isotherms (in unit erg/g). 

As a result, we may conclude that the contributions from pressure ionization and Coulomb 
coupling smooth the variation of the pressure and internal energy in the region of pressure ionization 
even when pressure ionization takes place in a very narrow density range. 
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5.2. Comparison of the EOS 

The EOS formulated in this paper has simple and explicit expressions for its practical use. 
However, its reliability requires comparison with the EOS data calculated by other EOS. Of many 
EOSs commonly used in astrophysics we choice the EOS calculated by Saumon et al. (1995) to 
compare our EOS. Their EOS is based on careful study of nonideal interaction and new physical 
treatment of partial dissociation and ionization by both pressure and temperature effects (see 
Saumon and Chabrier 1991,1992). They also made comparison with other EOS for a wide p — T 
domain, and gave detailed discussion and analysis. 

In Fig. 7 and Fig. 8 we plot the pressure and internal energy density on isotherms from 
our EOS and the Saumons' EOS for a pure hydrogen fluid. The difference at low density and high 
temperature comes from that our EOS includes the contribution from radiation while the Saumon's 
does not. To fit with Saumons' EOS on the internal energy density at very low temperatures, we 
have to ignore the zero temperature term of vibrational energy of molecular hydrogen, i.e., the term 
(l/2)/c6 v in equation (33). Comparison in the pressure in Fig. 7 shows difference mainly at low 
temperatures logT < 4.5 in the density range — 1 < logp < 1. Comparison in the internal energy 
density in Fig. 8 also shows the difference in the same region. These differences come from the 
different formalisms for the models of pressure ionization and Coulomb coupling. In spite of this, 
both of them prove the same tendency of the pressure and internal energy from neutral fluids to 
ionized plasmas. 



Fig. 7. — Comparison of pressure isotherms (in unit dyn/cm 2 ) between our EOS (solid lines) and the 
Saumons' EOS (circles) for a pure hydrogen fluid. The isotherms are (from bottom to top): logT = 
2.5, 2.9, 3.3, 3.7, 4.1, 4.5, 4.9, 5.3, 5.7. 

Another distinguished difference is revealed in the internal energy isotherm at temperature 
logT = 4.5 in Fig. 8. The difference covers a density range —4 < log/3 < — 1 and is specified by 
A log E = 0.07 or 17 percent. In this region our EOS results are more close to the EOS of MHD of 
which a difference is found with about the same amount (see Saumon et al. 1995). The difference is 
argued by Saumon et al. (1995) to come from a higher degree of ionization caused by the increase 
of the finite size of atoms with temperature and density dependence. 

For higher temperatures logT > 5, our EOS has almost the same result as Saumons' EOS in 
the internal energy. However, both EOSs have smaller values than the MHD EOS. This fact has 
been revealed by Saumon et al. 1995. This is because both EOSs do not take into account the 
Stark ionization while the MHD EOS does and yields higher degrees of ionization for a partially 
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Fig. 8. — Same as Fig. 7 but for internal energy density isotherms (in unit erg/g). 
ionized gas. 



6. Conclusions 

We present a single and complete algorithm of the EOS from cool-dense fluids (T > 10 2 K) to 
hot gases for a H-He mixture. The physical processes of pressure ionization and Coulomb coupling 
are simply modeled for its validity from cool-dense fluids to hot gases and for mixed elements. 
This allows to archive an EOS with explicit expressions which can be easily incorporated into the 
modeling of stars or planets. In addition, the simple and comprehensive algorithm also allows to 
incorporate more physics processes or add more elements into the EOS formalism. 

We show that the thermodynamic functions remain smooth even at the abrupt occurrence 
of pressure ionization at very low temperatures. The contributions from the Coulomb coupling 
and pressure ionization to the thermodynamic function compensate the effect of the rapid increase 
of free electrons when pressure ionization takes place. Although the quantum exchange effect of 
electrons is not included, the resulted algorithm produces comprehensive data of the EOS from hot 
gases to cool-dense fluids. 

This work was supported by the Science Foundation of Yunnan and the National Nature 
Science Foundation of China. Thanks are also due to D. Saumon and H. Chabrier for their EOS 
data which we accessed by internet ftp. 
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